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Abstract 

The convergence of the linear 5 expansion for the connected generating functional 
of the quantum anharmonic oscillator is proved. Using an order-dependent scaling 
for the variational parameter A, we show that the expansion converges to the exact 
result with an error proportional to exp(— a/V 1 / 3 ). 
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1 Introduction 



The linear delta expansion (LDE) is an analytic approach to field theory which has been ap- 
plied to a variety of problems (see for example Ref. |]T]). The approach is non-perturbative 
in the sense that a power series expansion is made in a parameter S artificially inserted into 
the action which interpolates between a soluble action Sq and the action for the desired 
theory S. The action is written 



S 5 = (l-S)S + SS, 



(1) 



where Sq contains some dependence on a variational parameter A. The generating func- 
tional for the theory may be evaluated as a power series in 5, which of course is usually 
evaluated to finite order. When 5 is set equal to one at the end of a truncated calculation 
of, say, a Green function G, some dependence on A remains, which is where the variational 
procedure makes its appearance. One such procedure is the principle of minimal sensitivity 
(PMS) H] which requires A to be a stationary point of the truncated Green function G^: 



dG N (X) 
dX 



0. 



(2) 



The PMS can provide both non-perturbative behaviour and convergence. It has been 
shown to ensure the convergence of the linear delta expansion for the vacuum generating 
functional for 4 theories in both zero || and one || dimension. For the finite-temperature 
partition function Z^{(5) at odd N there is only one stationary point, a global maximum, 
so there is no ambiguity. Problems can, however, arise when there are multiple solutions 
to the PMS condition Eq. (0), as is the case H for the connected vacuum generating 
functional W = InZ in zero dimensions. There, rather than applying the PMS directly 
to W, an order-dependent scaling for A was chosen which guarantees convergence for W. 
This was in fact the scaling resulting from the application of the PMS to the calculation 
of Z. In the present paper, we use a similar choice of scaling to provide a convergent 
series of approximants for the connected vacuum generating functional of the anharmonic 
oscillator. 

The simplest choice of delta expansion used in Ref. [|]] (using as trial action a free 
action with variable mass) gives the finite-temperature partition function as a function of 
z, the complex extension of 5: 



Z(z) = — / Dx exp 

Zq Jx(0)=x(J3) 



rP /l 1 

/ dr -x 2 + -(m 2 + 2gX)x 2 + zg{x A 
Jo \2 2 



Xx 2 



(3) 



where Z is the partition function for g = 0. It is this system that we consider here. The 
scaling required for the PMS which also gives convergence was shown in Ref. Q to be 
A oc N 2 / 3 . In the present paper, we use this scaling to prove the convergence of the 5 
expansion for W = In Z, provided that m 2 > 0. This scaling is not strictly PMS, simply 
a choice to give guaranteed convergence. 
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In the following Sec. |2], we review the proof of convergence in the zero dimensional 
model. This serves as a detailed preview of the arguments used in the main sections of 
the paper. In Sec. |3|, we set out the asymptotic evaluation of the remainder Rn{z) = 
Z(z) — Zn(z) for the anharmonic oscillator. The analysis is closely related to that of 
Ref. J3J. In Sec. f| we discuss the evaluation of Eq. ([$]) for Z(z) by saddle point methods 
and the cancellations that govern the convergence of Wjv- We compare our analysis with 
numerical calculations in Sec. |5|, and summarize our results in the conclusions, Sec. || 

2 Illustration of the method 

The method of the present paper is similar in spirit to that of Ref. ||, which considers the 
zero dimensional analogue of the vacuum generating functional, namely: 

roo 

Z®(z) = / dx exp(-[A(l - z)x 2 + gzx 4 ]) (4) 

J — oo 

An expansion of in powers of z was proved in Ref. |§ to converge to the correct value 
forZM(l) with an error like e cN . The PMS condition guarantees convergence and requires 
A to scale as \/N. 

In Ref. [|J it was shown that the convergence of Wn = {In Z n (z)}n\ z= i, where {/(^)}jv 
denotes the Taylor expansion of f(z) truncated at 0(z N ), is governed by the position of 
the roots of Z^(z). If any of these roots lies inside the circle \z\ < 1, the expansion for 
Wn(z) is not convergent, and in fact the position of the smallest modulus root determines 
the rate of convergence of the series for W^. The difficulty in such a proof is that while 
it is a relatively simple matter to calculate the coefficients c n in the series Z^(z) = T,c n z n 
(using simple integrals in zero dimensions or the methods of Caswell [f| and Killingbeck 
in one dimension), finding the roots analytically is highly non-trivial. 

In order to circumvent this difficulty, we write 

Z N {z) = Z{z) - R N {z) (5) 

From this point of view, zeros of Z^(z) arise from a cancellation between Z(z) and Rn(z), 
each of which can be evaluated by asymptotic methods. In Ref. ||, the remainder Rn{z) 
is bounded quite simply, while Z(z) is estimated by stationary phase integration. 

The saddle points of the exponent S in Eq. (4) are x = and x± = ±A(1 — z)/2gz. 
The mechanism identified in Ref. || for roots to lie inside the circle of mean radius \z\ = 
\cq/cn\ 1 I n was cancellation between the saddle-point contributions from x and x±. This 
can only occur along the Stokes line T defined by Re S(x±) = 0. However, part of T lies 
inside \z\ = 1, and cancellation along this part of the curve would have been disastrous 
for the proof of convergence. Such cancellations were ruled out in Ref. || by explicitly 
constructing the stationary phase paths for different values of z. It was shown that for 
\z\ < 1 the correct path did not pass through x±, in contrast to the region \z\ > 1, where 
all three saddle points were traversed. 
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Because this argument is difficult to generalize to a genuine path integral, we note here a 
simpler argument ruling out dominant contributions from x± in the region Re z < 1. In fact 
Z^(z) may be evaluated exactly as a modified Bessel function Ki((X(l — z)) 2 /(8g 2 z 2 )), 
provided that Re z < 1. The asymptotic behaviour of this function as A — > oo has 
no exponential dependence on iV in this region. Thus saddle points potentially giving 
a greater contribution than the Gaussian saddle point xq are excluded. The argument 
immediately generalizes to the left half plane without the need for contour rotation. The 
analysis in the present paper is similar, though we are unable to escape the need for contour 
rotation. 

In the region Re z > 1 all three saddle points can occur, with possible cancellations 
along the curve T. In this way it was shown that the smallest zero of occurs at 
Zm i n = a + \/a 2 — 1 where a = 1 + 3iri/(Na), with a = 1.3254.... This tends to unity 
sufficiently slowly as iV increases to ensure that W$ converges to with an error of 
the order of exp(— y^riV/a). 

3 Convergence of the remainder Rn{z) 

The remainder = Z(z) — Z^{z) can be written as 

If \ rP fl I \ 

R N (z) = — I Dx exp - / dr \-x 2 + -(m 2 + 2g\)x 2 + zg(x 4 - Xx 2 )) 

Z Jx(o)=x(/3) Jo \2 2 / 

x [1 - Q N (y)} , (6) 

where Qn(v) '■= e~ y {e y }N and y = —gz J dr(x 4 — Ax 2 ). 
By integrating the identity for 6 at given in Ref. |§: 

d ii N e~ v 

we can write 

l-Q N (y)= / dy' y 



lo * N\ 

Since z and hence y is in general complex, we must make a change of variable dependent 
on the sign of the real part of y. For Re y > , we may write y = \y\e l9 and change 
variables using y' = use %e 

r\v\ uj n 

1 - Q N (y) = J duo — exp (-ooe 10 + i(N + 1)0) . (9) 

For Re y < 0, we write y = \y\e i{e+7r) (— k/2 < 6 < vr/2) and use y 1 = ue i{e+n) : 

r\v\ uj n 

l-Q N (y) = ^ duj — exp(-uje i9 + t(N + l)(e + n)) (10) 
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Applying these to the remainder Eq. (|6]) for Re z > 0, we define weak and strong field 
regimes respectively by 



P 



P 



drx 



A I drx 2 > 
Jo JO 

rP rP 

A / drx 2 < / drx 4 
Jo Jo 



We divide the remainder Rn{z) into weak and strong field contributions and B^, and 
find 



A N (z) 



B N (z) 



— f Dx exp — f dr ( -x 2 H — (m 2 + 2g\)x 2 + g\z\ cos6* (x 4 — Ax 
Zq J a L J \2 2 

r\v\ u N 
x / da; — r exp (—a; cos 9) 
Jo N\ ' 

— I Dx exp — dr ( -x 2 H — (m 2 + 2^A)x 2 + g\z\ cos ^ (x 4 — Ax 2 ) ) 
Zq Jb L 7 \2 2 / 

/•Is/I a;^ 
x y aa; exp(a; cos 6 1 ) . 



12) 



(13) 



/?v( :)l < D\z\ N+1 exn (-ciV 2/3 ) . (14) 



The two regimes are dealt with separately in the following subsections. For Re z < 0, 
we simply interchange the role of the weak and strong field contributions to Rn{z). The 
calculations give a result similar to the zero dimensional case, namely that we may bound 
the remainder by 

\*In{z)\ < D\z\ N+1 exp 

The detailed derivation of Eq. (|14D for \z\ — 1 and its reliance on the PMS scaling A = 
{2^ 2 N 2 / g) 1 ^ with 7 = 0.186 are given in Ref. 0. For \z\ ^ 1 we see that Rn{z) is 
convergent for 

\z\<l + 0(N- 1/3 ) . (15) 

3.1 Strong fields 

This is the simpler regime to consider, as the analysis almost exactly mirrors that of Ref. M . 
Since the u integrand in Eq. (|l^) is monotonically increasing, we may simply bound by 
the value at the upper limit: 



B N (z)\ < 



1 



ZnNl 



f Dx exp — f dr f -x 2 H — (rn 2 + 2g\)x 2 + g\z\ cos 6 (x A — \x z 

J B J \ 2 2 



x l?/!^" 1 " 1 exp (\y\ cosi 



(16) 



Now using \y\ = g\z\ J dr [x 4 — Ax 2 ], we find that 



B N (z)\ < 



Z a N\ 



x 



- / Dx exp — dr l -x 2 H — (m 2 + 2gX)x" 
. J b «/ \ 2 2 

9 J dr (x 4 — \a 



(17) 



which is just l^l^ 4 " 1 times the given in Ref. [|[]. Thus with no further analysis we are 
able to bound the strong field remainder by 

\B N {z)\ < \z\ N + 1 C(m,g){3N 4 ^ex P (-NS B ( 1 )) (18) 

with Sb and 7 as defined in Ref. J§]. 



3.2 Weak fields 

The weak field contribution to the remainder is given by 
/ Dx ej I . 1 . > 

'A 

'\v\ . u 



An{z)\ — — I Dx exp — / dr ( -x 2 H — (m 2 + 2g\)x 2 + g\z\ cos #(x 4 — Ax z 
^0 -A L j V 2 2 

rlvl 

x y duo -j^ exp (— tucos^) . 



(19) 



If we make a further change of variable to = g\z\a J dr[Xx 2 — x 4 ], the familiar factor l^l^" 1 " 1 
already emerges: 



IJV+1 ,1 



A N (z) 



Z n N\ 



- f da f Dx exp - / dr (-x 2 + -(m 2 + 2g\)xA 
! Jo J a 1 J \2 2 / 



x (g J dr (\x 2 — x 4 )^j exp cos 0(1 — a) J dr (Ax 2 — x 



(20) 



For the allowed range of 9, it is then clear that 

\N+1 ,1 



A N (z)\ < 



\z\ 

X 



— — - da Dx exp - dr [ -x 2 + -(m 2 + 2g\)x' 
oiV! jo ii L J \2 2 

/ /■ \ iV + 1 

/ L / dr(Ax 2 -x 4 )J e: 



g|z|(l — cr) / dr (\x 2 — x 4 ) 



(21) 



We may then further bound An(z) by removing the term (1 — a)gzx 4 from the last 
exponent in Eq. ([H]). We also use Stirling's approximation for the factorial. We follow 
Ref. [|J and allow the calculation to proceed for either sign of m 2 by writing: 



A N \ < 



where 



IAT+1 -1 



f Q daJ A Dx{gJdr(\x 2 -x 4 ))^ - J q dr Qx 2 + \ 



^2-kNZq 
x exp (—NSa(x, cr, z)) 



1 2 1 2 
\m x 



(22) 



Sa{x, a, z) 



N 



dr I C (<?A (1 - |z|(l - a)) + m 2 fl(-m 2 )) 



In f — f dr (Xx 
\N Jo V 



2 4 

x 



(23) 



which for z = 1 is equal to the Sa defined in Ref. [|J. 

The calculation proceeds in exactly the same manner as in Ref. |3J , with the definition 
Jo x 2 dr = (3XU and the use of the Cauchy-Schwarz identity to show that 



S A > all - \n(aU(l - U)) - 1 
with a and a the z ^ 1 extensions of those defined in Ref. [ 

g(3\ 2 a 

m 2 p\ 



a. 



N 



a 



ft \ \z + 



+ 



a 



N 



6{-m z ). 



We minimize the bound on Sa with respect to U and take the minimum at 



[7=1 + 1 
2 a 



which gives a similar result to that of Ref. [|J : 



Sa > F(a) - In ( - 



a 



where 



a. 



a 2 



F(a) = --Wl + T + ln|l 



It is true for the present definition of a that for A oc N 2 ^ 3 and large N, 



d_ 

da 



F(a) - In ( - 



< 



(24) 



(25) 



(26) 



(27) 



(2* 



(29) 



so that the integral over o in Eq. (|22]) is dominated by its value at the upper limit. Thus: 



liV+l 



AatI < 



V2irNZ J A 
x exp 



[g dr 



Xx — x 



exp 



• 1 1 2 2 

-X H 771 X 



N ( F(ft ) - In ( -ft 



with the subscript on ft and a denoting a = 1, giving 

^A 2 



(30) 



ft 

o: 



N 

m 2 f3\ 2 , 
ft + —jy—9(-m / 



(31) 
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Now if we use the relation: g Jq dr [Xx 2 — x 4 } < g(3X 2 /A and the fact that the functional 
integral in the numerator of Eq. ( |30|) is bounded above by that in the denominator, we 
obtain 

g(3X 2 \z\ N+l 
A N < ^ L— 
~ 4 v / 2^iV 

Using the scaling A oc N 2 / 3 we have a « a [l + 0(N~ 2 ^ 3 )]. For large N, a oc iV 1//3 and 
we may finally write 



exp 



—N ( F(a ) — In 



2 J 



(32) 



JV+1 



AT5/6 



exp 



(33) 



4 Saddle point expansion of Z(z) 

The partition function may be written as 

Z{z) = ^( Dxe*v(-S[x}) (34) 

where 

S'H := f dr (-x 2 + -M 2 x 2 + gzx 4 ) (35) 
Jo \2 2 / 

and M 2 = m 2 + 2g\(l — z), which is large for A oc N 2 ^ 3 and z ^ 1. 

Performing a saddle point expansion of the path integral requires finding the solutions 
of the Euclidean equation of motion 

x = M 2 x + Agzx 3 , (36) 

with x(0) = x(/3), which has the static solutions: 

x = Xq =0 (37) 

x= x - =± t W* (38) 

There is also a non-static instanton solution. This is discussed by Zinn- Justin H in the 
case of negative g and real M 2 > 0. This class of solution also arises in the strong coupling 
contribution to the remainder in Ref . || , the modified action giving an equation of motion 
equivalent to negative coupling. This solution is 



—M 2 1 

Xl{T) = V^7 cosh(M(r-r )) ' (39) 

where To must strictly be taken as r = |/3 to satisfy the periodic boundary conditions. 
We note that this solution takes the form 1/ cos(M(r — r )) for real z, such that z > 
1 + m 2 /(2gX). For general complex z the solution is an hybrid. 
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A more general class of solutions is discussed by Richard and Rouet for the double 
well oscillator. The action of Eq. (|3~5j) may be rescaled to the form of Ref. which is: 



Srr[x] = [ dt 

J-T 



T 

(40) 



• 2 , 1/2 1 \2 

~x H — x — 1 



.2 2 

The solutions to the modified equation of motion take the form 

x RR (t) =a(l + — — Q 7 3 '-j — ] (41) 

where a is an integration constant and V(t + u\u, u>') is a Weierstrass elliptic function flTO] 



whose periods a;, a/ are determined by the boundary conditions and form of the elliptic 
integral inverted by Eq. (f|l]). The cosine instanton solution discussed above is equivalent 
to Eq. ( f4"T|) with uj'/uj — > zoo. 

The classical action for all the saddle points at large M may be evaluated simply: 

S[x ] = (42) 
M 4 8 

S[x ± ] = -— ^ (43) 
logz 

M 3 

= -- (44) 
3[xrr] = S[x±] + kS[x T } , (45) 

where k is a positive integer, except that in the case u'/u — > zoo, cancellations occur in 
the action of Richard and Rouet to give S^x^r] = S^x/]. 

Having evaluated the classical actions, we must evaluate the regions of the z plane in 
which the various saddle points are dominant. These regions are illustrated in Fig. 1 for 
the numerical values m — 1, g — 1/2, [3 = 2 and A = 11.619 (N = 75). The solid curve 
T is the curve Re S[x±] = and the dotted curve Tj is the curve Re S[xj] = 0, while the 
dashed vertical line marks Re z = 1 + m 2 /2g\. All three meet at the point z = 1 + m 2 /2gX. 

Inside Tj, the Richard and Rouet saddle points are dominant. Inside the loop of T, the 
saddle point x± dominates over xo, though this is unimportant given the overall dominance 
of S[xrr}. Outside T/ the zero saddle point becomes dominant over all its competitors. 
This remains the case up to the boundary T to region C, in which x± is dominant. Neither 
Xj nor xrr can dominate the integral in this region. 

It would then appear that instantons dominate for small \z\, Xq is dominant for in- 
termediate | z |, and for large \z\ the saddle points x± are the important ones. However, 
as we learnt from the zero- dimensional case, it is important to check whether apparently 
dominant saddle points do in fact contribute. In zero dimensions Ref. || it was possible 
to trace out the stationary-phase paths and to see that these did not pass through the 
non-zero saddle points when \z\ < 1. For a path integral, such a procedure is extremely ill 
defined. However, we are able to use instead an analogue of the Bessel function analysis 
given above for the zero dimensional case by realizing that for Re z < 1 + m 2 /(2g\) the 
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J Dx exp 


rl! dT ( 









action is generically a real single well, the complex part of the mass term simply adding a 
phase. Provided then that Re(gzx 4 ) > the partition function is bounded above by that 
for the pure harmonic oscillator. The harmonic oscillator partition function is given by 

+ Ir 6 (mV 

= C(det(Sf + ReM*))~5 , (46) 

where C is a normalization constant. The Gaussian integration about the saddle point 
x = in Eq. (f51) gives: 

\Z\ = C\det(d? + M 2 )\-^ . (47) 

We are thus able to state that in the region < Re z < 1 + m 2 /(2g\) there is in fact no 
exponentially increasing contribution from saddle points with Re S[x] < 0. The zero saddle 



point is the only one which can contribute. From the form of Eq. fll4), it is then clear that 
no zeros of Zn(z) are possible for \z\ < 1 + 0(A^ -1 / 3 ) in the region < Re z < 1+m 2 /(2g\). 
In fact the zeros will be expected to occur on a ring of radius l + O^- 1 / 3 ). 

The left half of the z plane may be included by rotating the contour of integration to 
make R,e(zx 4 ) > while keeping both Re(x 2 ) > and Re(M 2 x 2 ) > 0. This procedure is 
similar to that discussed in Ref. ||. In the upper half plane, we rotate the contour by an 
angle = (jr/2 — 9)/4, where 9 = argz. In the lower half plane, we simply rotate the 
contour by an equivalent amount, but in the opposite sense. 

In the region Re z > 1 + m 2 / (2g\), we need only consider the stationary saddle points 
xq and x±. In the region B, xq dominates and the requirements for cancellation are the 
same as those discussed above. No zeros affecting the rate of convergence of Wn can arise 
in this region for m 2 > 0, and we expect the ring of zeros to continue smoothly into this 
region. For the non-Borel summable case m 2 < 0, this proof breaks down, since zeros of 
Zn(z) could now occur for \z\ < 1. 

In region C, where x± is the dominant saddle point, Z(z) contains exponential depen- 
dence on N, allowing cancellation with R^(z), which grows with |*| for |*| > 1 + 0(NV 3 ). 
The zeros in region C will thus tend to be further out than those in B. 

An exception to this is the region close to T, the Stokes line where Re S[x±] = 0. 
There cancellations can occur between the two saddle-point contributions from xq and x±, 
giving a smaller effective exponent for Z(z), and the possibility of roots of smaller modulus. 
However, since this branch of T lies to the right of Re z = 1 + m 2 / (2g\), such zeros again 
do not affect the convergence of Wn for m 2 > 0. Having established that all the roots 
satisfy \z — 1| > m 2 /2g\ = 0(l/iV 2 / 3 ), the analysis of Ref. || shows us that W N tends to 
the correct value W with an error 

n N := W-W N 

oc exp (-ciV 1/3 ) (48) 
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5 Comparison with numerical results 



The above analysis of the position of the roots of Z^{z) can be checked numerically. 
Figs. 2, 3 and 4 show the roots of Zn(z) for (3 = 2, m = 1, g = 1/2 and N = 25, 45 and 
75 respectively. In each case, A is chosen as the PMS value for Z^(z). In the previous 
section, it was noted that the smallest root occurs where Re S[x±] = 0. In order to make 
the comparison, we have plotted the curve T on each graph. 

In Fig. 2 we see that at order 25, the roots lie to a good approximation on a ring, 
the radius of which is determined by the ratio of the first and last coefficients in the 5 
series. Fig. 3, for order 45, shows that two pairs of roots have broken away from this ring, 
one pair moving out and one pair moving in. The inner pair of roots have the smallest 
modulus, and lie sufficiently close to T to indicate that the mechanism outlined above is 
indeed operative. At order 75 in Fig. 4, two pairs of roots are tracking in along the curve 
T. This behaviour very closely follows that described in Ref. || for the zero dimensional 
case. 

We have checked the convergence of the expansion for W numerically. Taking A at as 
the unique PMS value for Z^ we obtain Table 1, showing convergence to 15 significant 
figures by N = 85 for both (3 = 2 and (3 = 5. The convergence is not monotonic, as was 
the case for Z N ; rather it oscillates around the exact value with exponentially decreasing 
amplitude. 

Notwithstanding the problem of multiple PMS points for Wn(X) itself, it is interesting 
to explore the dependence on A of Wjq. For the values of N we have considered there seems 
to be an initial, extremely broad maximum, which gives a very accurate estimate of W, 
followed by a series of secondary maxima and minima of progressively decreasing accuracy. 
For N = 45 and (3 = 2 the value at this first maximum is -1.353868180362611, differing 
from the exact value only in the 16th. significant figure. A similar situation occurs for 
(3 = 5. 



N 


W N ((3 


= 2) 




W N {(3 = 5) 


17 


-1.353* 


36* 


32048367371 


-3.4* 


305935670797027 


25 


-1.353* 


m 


31* 


W4666174 


-3.4* 


305880725830912 


35 


-1.353* 


m 


m 


W3659796 


-3.4* 


305878403631746 


45 


-1.353* 


m 


m 


503626513 


-3.4* 


305878324251261 


55 


-1.353* 


m 


m 


503625918 


-3.4* 


305878320234629 


65 


-1.353* 


m 


m 


503626148 


-3.4* 


305878319968700 


75 


-1.353* 


36* 


31* 


503626174 


-3.4* 


305878319947731 


85 


-1.353* 


36* 


m 


503626168 


-3.4* 


305878319945885 


Exact value 


-1.353* 


36* 




503626170 


-3.4* 


305878319945603 



Table 1: Results showing convergence of 5 expansion for Wn 
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6 Conclusions 



In this paper, we have proved the convergence of the optimized 5 expansion for W = In Z, 
where Z is the finite-temperature partition function of the anharmonic oscillator. With 
the variational parameter A chosen to scale with the order iV of the expansion as 

A = (2 1 2 N 2 /g) 1/3 (49) 

we found that 

W - W N oc exp (-ciV 1/3 ) . (50) 

Apart from the fact that in a field theoretic context it is generally the connected Green 
functions which are the relevant quantities, the original motivation H for looking at W 
was the non- uniform convergence with /3 of the delta expansion for Z^, in particular in the 
limit f3 — > oo. Since this non-uniformity corresponds to the limit of large spacetime volume, 
it seemed plausible that the difficulty would be much less severe for W, the generating 
function of connected diagrams, which should depend only linearly on the (large) volume 
cut-off. However, the present proof of convergence for Wn does not address this problem, 
since it relies as an intermediate step on the convergence of the sequence Zn- It remains 
an open question whether the PMS criterion applied directly to Wn itself, in spite of the 
problems of multiple stationary points, can give rise to uniform convergence as (3 — » oo. 

In the meantime a paper by Guida, Konishi and Suzuki [II] has appeared, which, using 



a completely different approach involving a dispersion relation in g, has proved convergence 
of the 5 expansion for the individual energy levels of the anharmonic oscillator for a variety 
of scaling exponents rj in A^r oc N* 1 . The PMS choice we have used here, with rj = 2/3, lies 
at the edge of their range 2/3 < r] < 1. It is possible that their method could be extended 
to a proof of convergence for W, or the free energy F = W//3. 

The proof of Ref. |L1| is restricted to m 2 > 0, as is the case here, even though it was 
shown in Ref. [|J that for Z the sequence of approximants Zn converges for either sign 
of m 2 . As a general principle the field theory about which we expand should capture 



as closely as possible the essential features of the system under investigation | I2]] . For 
the double-well oscillator it may well be that to obtain a convergent expansion for the 
connected Green functions a more sophisticated trial action is needed than the one used 
here and in Ref. [II], namely a free action with positive m 2 . 

The extension of the method to higher dimensions should be possible, since only saddle 
point techniques are used. In higher dimensional field theories we must, however, take 
account of the interplay between the 6 expansion and the renormalization procedure. Such 
a scheme has been successfully applied to the Gross-Neveu model in the large N limit 
using the 8 expansion [13| and a related scheme [14 and 4 theory in four dimensions in 



the Gaussian approximation [[US]]. 
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